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Abstract 

Classical acoustoelasticity couples small-amplitude elastic wave propagation to 
an infinitesimal pre-deformation, in order to reveal and evaluate non-destructively 
third-order elasticity constants. Here, we see that acoustoelasticity can be also be 
used to determine fourth-order constants, simply by coupling a small-amplitude 
wave with a small-but-finite pre-deformation. We present results for compressible 
weakly nonlinear elasticity, we make a link with the historical results of Bridgman 
on the physics of high pressures, and we show how to determine "D" , the so-called 
fourth-order elasticity constant of soft (incompressible, isotropic) solids by using 
infinitesimal waves. 
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1 Introduction 



Acoustoelasticity is now a well established experimental technique used for the non- 
destructive measurement of third-order elasticity (TOE) constants of solids, and its prin- 



ciples can be found in standard handbooks of physical acoustics, such as Pao et al. (1984); 



Kim and Sachse (2001). The underlying theory, however, is quite intricate, and over the 
years several authors have produced different and irreconcilable expressions for the shift 
experienced by the wave speed when elastic wave propagation is coupled to a pre-strain. 
A common mistake found in the literature consists of the implicit assumption that since 
a small pre-strain and a small-amplitude wave are described by linearized equations, they 
can be superposed linearly. With that point of view, the coupling between the two phe- 
nomena is of higher order and can be neglected in the first approximation. The flaw in 
that reasoning is simply that these phenomena are successive, not linearly superposed. 
Experiencing an infinitesimal pre-strain from a stress- free configuration is a linear process; 
propagating an infinitesimal wave in a stress-free configuration is another linear process. 
But in acoustoelasticity, the wave travels in a pre-stressed, not a stress-free solid, and 
the laws of linear elastodynamics must be adapted to reflect this fact. The final outcome 
of this analysis is that not only second-order, but also TOE constants appear in the 
expression for the speed of an acoustoelastic wave. 

For example, a longitudinal wave travelling in a solid subject to a hydrostatic stress 
al will propagate with the speed Vql given by 



7A + 10/i + 2A + lOfi + 6C 
3A + 2/i 



XI] 



where p is the mass density in the unstressed configuration, A and fi are the (second- 
order) Lame coefficients, and A, B, C are the (third-order) Landau coefficients ( Landau] 
and Lifshitz, 1986). This equation was correctly established as early as 1925 by Brillouin 



and confirmed in 1953 by Hughes and Kelly, although several erroneous expressions have 



appeared in between and since (see, e.g.. Birch (1938); Tang (1967)) 



From our contemporary perspective, the easiest way to re-establish and further this 
expression is to rely on the modern theory of incremental (also known as small-on-large) 
elasticity, which presents compact expressions for the tensor of instantaneous elastic mod- 
uli, A{), given by its components 



(A 



Oiijj 



A-Oijij 


= J 


A-Oijji 


= J 


A-Oijij 


= J 


, Aoijji 


= J 



A? -A? 



i ^ j, Xi ^ A 



J' 



_,XjWi-XiW, 



i ^ j, Xi ^ X 



[1.2) 



Xj, 
Xj, 



in the coordinate system aligned with the principal axes of pre-strain, which, since the 
material is isotropic, coincide with the principal axes of pre-stress. Here, W is the strain- 
energy density per unit volume, which is a symmetric function of the principal stretches 
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Ai,A2,A3 of the deformation, J = A1A2A3, Wi = dW/dXi and Wij = d'^W/dXidXj, i,j G 
{1,2,3}; see, for example, Ogden (1997 2007[ ) for details. The corresponding principal 
Cauchy stresses are given by 0", = J~^XiWi, i = 1, 2, 3. 



It is worth noting here for later reference that it follows from (1.2)2 that 
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A small-amplitude body wave may travel at speed v in the direction of the unit vector 
n with polarization in the direction of a unit vector m provided the eigenvalue problem 



Q{n)m = pv'^m 

is solved with pv"^ > 0, where Q{ti) is the acoustic tensor, which has components 



Qij{n) = J A 



[lA) 



[1.5) 



Note that sometimes the acoustic tensor is defined by (1.5) without the factor J, in which 



case the p in (1.4) would be replaced by the density in the deformed configuration, i.e. 

j-V- 

Notice that in this theory the amplitude of the acoustic wave is infinitesimal, but 
that there is no restriction on the choice of W or on the magnitude of the underlying 
pre-strain because we are in the context of finite nonlinear elasticity. It follows that the 
expressions can be specialized in various ways, in particular to weakly nonlinear elasticity 
(where W is expanded in terms of some measure of strain) and to special pre-strains. 

Hence, to access the TOE constants, we take W as 
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where /„ 



tr(^;" 



n 



1, 2, 3, are invariants of the Green-Lagrange strain tensor E. 
The components (1.2) can then be expanded to the first order in terms of, for example, 
the volume change e = J — 1 m the case of a hydrostatic pressure, or of the elongation 
d = Ai — 1 in the case of uniaxial tension. 



To access fourth-order elasticity (FOE) constants, we take W as 
W = WpoE = Wtoe + EIJ3 + Fllh + Gil + Hll 



[l.l) 



where E, F, G, H are the FOE constants and we push the expansions of the wave speed 
up to the next order in the strain. This is what we refer to as the large acoustoelastic 
effect. In fact, the main purpose of this investigation is to provide a theoretical backdrop 
to the current drive to determine experimentally the FOE constants of soft solids such 
as isotropic tissues and gels in order to improve acoustic imaging resolution (see, e.g.. 



Hamilton et al.| (|2004D; [Zabolotskaya et al.| (|2004[ |2007D ; [Gennisson et al. (|2007D ; Renier 
et al.| ( |2007D ; [Jacob et al.| ( |2007D ; [Renier et al.| ( |2008b|aD ; [MiTonov et al.| ( |2009f ) 



Because soft solids are often treated as incompressible, so that the constraint J = 1 
must be satisfied, we shall also consider the expressions for the TOE and FOE strain- 
energy densities in their incompressible specializations, specifically 



W 



(1.8) 
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for TOE incompressibility and 



W = iih + -h 



:i.9) 



for FOE incompressibility. Hence, in the transition from compressible to incompressible 
elasticity, the number of second-order constants goes from two to one, of third-order 
constants from three to one, and of fourth-order constants from four to one. This was 
established by Hamilton et al. in 2004, although it can be traced back to Ogden in 
1974 and even earlier, to Bland in 1969. Specifically, in that transition we note that the 



elasticity constants behave as (Destrade and Ogden, 2010) 



A — J- OO, — )■ -Ey/S, 

for the second-order constants, where Ey is the (finite) Young's modulus, 

A//i = 0(l), i?//x = 0(A/^), C/^i = 0{\^/^?) 
for the third-order constants. 



E//i = 0(A//i), F//x = 0(AV/i^), G//i = 0(A//x), H/fi = 0{Xyfi'' 
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for the fourth-order constants, and {C + F)/n = 0{X/n) for a combination of second- 
and third-order constants. Some constants have the explicit limiting behaviour 



5/A^-l, E/\^A/3, G/A^l/2, B + G + \/2^ D, 
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where D remains a finite quantity, of the same order of magnitude as fi; see Destrade] 
and Ogden (2010) for analysis of the behaviour of A, B, . . . , H in the incompressible 



limit and the connection with D. Note that the limiting behaviour for E in terms of the 
initial Poisson's ratio and Young's modulus Ey is thus (1 — 2h')E — )■ 4:Ey/9. This was 



shown in Destrade and Ogden (2010), but mistyped in equations (49) and (89) therein 
as (1 - 2iy)E AEy/3. 

In this paper, we treat in turn the case of hydrostatic pre-stress (Section |2]) and of 
uniaxial pre-stress (Section [3]), and we provide expansions of the body wave speeds up to 
the second order in the pre-strain, and also in the pre-stress, for compressible solids and 
in the relevant incompressible limits. 



2 Hydrostatic pressure 

Consider first a cuboidal sample of a compressible solid with sides of lengths (Li, L2, L^) 
in its (unstressed) reference configuration. We define the reference geometry in terms 
of Cartesian coordinates (Xi,X2,X3) by < Xj < Lj, i = 1,2,3. The material is then 
subject to a pure homogeneous strain xi = AiXi,a;2 = X2X2,X3 = A3X3 and deformed 
into the cuboid < Xi < k, i = 1,2,3, where {xi,X2,X2,) are the Cartesian coordinates 
in the deformed configuration, and the constants Ai,A2,A3 are the principal stretches 
of the deformation. We now specialize the deformation to a pure dilatation so that 
Al = A2 = A3 = J^/^. Since the material is isotropic the Cauchy stress cr is spherical, 
a = al say, with a > (< 0) corresponding to hydrostatic tension (pressure). 
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The pre-stress is computed as a = J^^XiWi or, more conveniently, as a = W'{J), 
where W{J) = W{J^^^, J^^^, J^^^). Up to the second order in the volume change e = J— 1, 
we find 



where k = W"{1) = A + 2yu/3 is the infinitesimal bulk modulus, and 



(2.1) 



W"'{1) = -K+ -{A + 9B + 9C). 
9 



(2.2) 



(It is a simple matter to check that the expansion (2.1 ) is equivalent to one first established 
by Murnaghan (1951).) Conversely, the volume change is expressed in terms of the 
hydrostatic stress as 

' ' (2.3) 
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-a 



K 2^3 

Here the pre-stress does not generate preferred directions and the solid remains 
isotropic. It follows that two waves may propagate in any direction, one longitudinal, 
with speed f^, and one transverse, with speed fy, given by 



pvl = J Ann, pv^ = JA 



■01212; 



(2.4) 



where these quantities are computed from the formulas (1.2). Expanding in terms of the 
volume change, we obtain 



where 



pvj^ = A + 2/i + glE + b^e'^, pv^ = p + axe + bxe^ 



flL = ^A + y/i+ y5 + 2C, aT = \ + 2p+^A + B, 



(2.5) 
(2.6) 



are the coefficients for the classical (linear) acoustoelastic effect, and 

18 9^^ 9 9 3 3 3 9 

h^ = -\+-^+-A+—B + C + E + F + -G (2.7) 
2 2 6 3 ^ ' 

are the coefficients of the large (quadratic) acoustoelastic effect. 

Alternatively, we may use (2.3) to express the wave speeds in terms of the pre-stress 
rather than the pre-strain, as 



pv\ = X + 2p + clct + dicr'^, pv^ = p + CtC + drcr'^, 
thus recovering the classical formulas 



cl = — , Ct = — , 
K K 



(2.8) 



(2.9) 



for the (linear) acoustoelastic effect (Hughes and Kelly, 1953), and establishing the for- 
mulas 
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Figure 1: Fit of the data measured by Hughes and Kelly (1953) for pyrex subjected to 
hydrostatic pressure P (in bars) based on the equations (2.8) and (2.11) with a = —P 
describing the large acoustoelastic effect: M and K are expressed in units of 10^ bars. 



for the large (quadratic) acoustoelastic effect, where W"'{1) is given in (2.2). 

For an application of these results, we turn to the classical data of Hughes and Kelly] 
(1953). They performed acoustoelastic experiments and plotted the stress-dependent 



shear and bulk moduli, defined as 



M 



PVt, 



K 



(2.11) 



against the hydrostatic pressure P = —a for polysterene and for pyrex. For the former, 
the variations are clearly linear, and nothing would be gained by including FOE effects. 
For the latter, there is a marked departure from the linear acoustoelastic effect at high 



pressures, and we thus use (2.8) to obtain a better fit to the data. By digitizing the data 



and performing a standard least square optimization, we found that for pyrex, 

M(P) = 2.73 - 8.45 x lO^^P - 1.21 x lO'^P^ 
K{P) = 3.27 - 7.80 X lO^^P + 4.07 x lO'^P^ 



(2.12) 



where M and K are expressed in units of 10^ bars and P in bars. The fit to the data is 
shown in Fig. [1} 



According to (2.8), the quantities (2.11) give access to A and /i, to two linear combi- 



nations of the third-order constants, and to two linear combinations of the fourth-order 



constants. Specifically, we may solve (2.11) to find the Lame constants as 

X = K{0)-^M{0), /i = M(0). 



(2.13) 
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Similarly, we find two linearly independent combinations of the three third-order con- 
stants: 



A + 3B = -3K{0) [1 + M'(0)] -4M(0), A + 9B + 9C = -9K{0) [1 + K'{0)] /2. (2.14) 
Finally, we also find two combinations of the fourth-order constants: 



E + F + lG=lKiO) + ^M{0) 

+ Im 



1 3 

K'{0) + ^M'(O) + K'{0)M'{0) + K{0)M"{0) 



E + 3F + G + 9H = ^K{0) + ^7^(0) 



K'iO) + I {K'{0)f + lmK"{0) 



(2.15) 



Hence, for the pyrex data: 

A = 1.45 X 10^ bars, 

A + 3B = -1.24 X 10^ bars, 

E + F + 2G/3 = -1.30 X 10^ bars. 



/i = 2.73 X 10^ bars, 

A + 9B + 9G = 1.00 X 10^ bars, 

E + 3F + G + 9H = 3.65 x 10^ bars. 



(2.16) 



Clearly, the data of Hughes and Kelly (1953) are not sufficient to determine all the third- 



and fourth-order constants. An additional relation is needed to determine experimentally 
all the third-order constants, and a further two relations to determine the fourth-order 
constants. These may be provided from, for example, wave speed measurements in uni- 
axially deformed samples, as we show in the following section. 

We conclude this section by evoking the work of P.W. Bridgman on "The Physics 



of High Pressures" (Bridgman 1945), which won him the 1946 Nobel Prize. He car- 



ried out countless high pressure measurements on solids, and fitted his data with a 
quadratic equation for the volume change, in the form AV/Vq = —aP + bP^. Direct 
comparison of this experimental law with (2.3) reveals the identifications a = 1/k, and 
b = —W"'{1)/{2k,'^). Further, acoustoelastic measurements give direct access to these 
constants, because k = K{0) and W'"{1) = —K{0)[2 + K'{0)]. In Bridgman's experi- 
ments, very high pressures are required to obtain second-order deformations; however, his 
coefficients a and b can also be determined by pressurizing the sample only linearly, and 
then measuring the speeds of transverse and longitudinal waves. In fact, had Bridgman 
been able to propagate such waves in his finitely deformed samples he would have had 
access to the third term in the expansion, AV/Vq = —aP + bP^ — cP^ say, because the 



continuation of (2.3) is 



AV 

^0 



--P-^,W"\1)P'--^^ 



3 (W"'{l)y - kW^''\1] 



(2.17) 



where iy*^'*''(l) is found to be expressible as 
iy(4)(i) 



19 4 8 

—K - -{A + 9B + 9G) + -{E + 3F + G + 9H), 



(2.18) 



in which each term can be determined experimentally using the equations of large acous- 
toelasticity above. 
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3 Uniaxial tension 



We again consider the cuboidal sample as in Section [2} but now the cuboid is subject to 
a uniaxial tension in the Xi direction so that it deforms homogeneously with elongation 
Ci = /i/Li — 1 = Ai — 1 in the xi-direction. By symmetry it contracts laterally and 
equibiaxially with elongation 63 = 62 = h/L2 — 1 = A2 — 1. The nominal stress (axial 
force per unit reference area) is Jcxi/Ai = aiXl- 



The condition that the lateral faces are free of traction is 



0^3 



(T2 = J-''X2W2 = 0, 



(3.1: 



which determines the extent of the lateral contraction — 1 < 62 < (assuming that 
A2 < 1) and can be written in terms of ei, up to second order, as 



where 
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X^A lX{X-2fi) 
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KA + /i)2 8(A + /i)3 L2(A + /i)2 j2(A + /i) 2(A + /i)3' 

(3.3) 

and again k = A + 2/i/3 is the infinitesimal bulk modulus. Note that ( |3.2 ) collapses to 
62 = — ei/2 + 36^/8 in the incompressible limits (1.10) and (1.12), as expected from the 



expansion to second order of the connection A2 = Aj^ 
incompressible solid in uniaxial tension. 



1/2 



between the stretch ratios of an 



Substituting (3.2) with (3.3) into the expression ai = J ^XiWi for the uniaxial stress, 

stress 



we obtain the relation between the pre-stress and the pre-strain, up to the second order, 
as 



X + fi 



ei + 7ei, 



where 

3nfi (5A + 3/i) 
2{X + fif 



A^ 



4(A + /i)^J 



^^3M3A^+4A^^+V)5 + 
2{X + fif 



(A + /i) 



(3.4) 



rC. (3.5) 



(Note that although 3K/i/(A + /i) = Ey, the infinitesimal Young's modulus, we shall not 
use Ey hereon.) Conversely, 



ei 



A + /i 



-ai 



A + /i 
3K/i 



(3.6) 



In the incompressible limits (1.10) and (1.12), the stress-strain and strain-stress relations 
reduce to 

^1 = 3/ie, + 3 (/i + ^) el e, = ^^1 - ^ (/^ + f) (3-7) 

Now we examine the possibility of a small-amplitude body wave travelling in the 
direction of the unit vector n with polarization in the direction of a unit vector m. For 
a general direction of propagation, the speeds of the different waves are found as the 



eigenvalues of the acoustical tensor Q{n) in (1.5), i.e. as the roots of a cubic. In this 
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paper, however, we focus on the wave speeds that are relatively simple to determine 
experimentally. The speeds of non-principal body waves are not easily accessible because 
they would require transducers to be placed at an angle to the faces of the cuboid, and 
not in full flat contact, which would lead to additional transmission problems. This leaves 
the principal body waves, as explained by Hughes and Kelly (1953): first, waves in the 
direction xi of the tension, and, second, principal waves in the X2 and X3 directions, which 
are equivalent by symmetry. 

In the direction of tension there exists a longitudinal wave with speed vu, say, and two 
transverse waves, propagating with the same speed vu, say (in fact, these latter two waves 
may be combined to form a transverse circularly-polarized wave). Hence, n = m = ei 
for the longitudinal wave, and n = ei, m = 62 for the transverse wave, where Ci and 
62 are unit basis vectors corresponding to the coordinates Xi and X2, respectively. Using 
(1.4) and (1.5) we then have 



^"^01111; P'Vi2 — </"4oi212- 

Expanding the elastic moduli in terms of the elongation ei, we obtain 

A + 2yU + aiiCi + buel, pvl^ = p + 01261 + &i2ei. 



(3.8) 



(3.9) 



where 



ail 



ai2 



{X + 2B + 2C) + 2 (2A + 5/i + A + 2B) 



for the classical acoustoelastic effect, and 
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, 9(A + 2/i) ^ , 
'^^ = '^ + ^(AT^^ + 
5A2 + 2A/i + 2/^2 



9/i 



-B 
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(3.12) 
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for the large acoustoelastic effect. 

For propagation perpendicular to the direction of tension, we take n = 62- Then there 
exists a longitudinal wave, with m = 62 and speed ^22, a transverse wave polarized in the 
direction of tension m = Ci with speed V21, and a transverse wave polarized perpendicular 
to the direction of tension, with m — and speed V23. These speeds are given by 



PV22 — •/^02222, P^il — •^^02121, PV23 — -'^^02323 ■ 

Expanding the elastic moduli in terms of the elongation ei, we obtain 

= A + 2/x + 02261 + 62261, pvjk = fJ' + a2kei + b2kel, k = 1,3, 

where 

2A(A + + A^ + 2(A - n)B - 2fiC 



0,22 — 
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A + // 

(A + 2/x)(4/x + A) +4/x^ 
4(A + A^) 
A(4/x + >l) - 2ixB 



2(A + //) 

are the classical acoustoelastic coefficients, and 
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(3.15) 
(3.16) 
(3.17) 
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(3.18) 



are the 'large acoustoelasticity' coefficients. 

We may also find expressions for the acoustoelastic effect in terms of the pre-stress 
(Ji as 

pvl = X + 2ij, + Ciiai + diial, pvf^ = IJ' + Cikai + dikaf, (3.19) 
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where i = 1,2 (no sum on repeated i), and k = 1,3 {k i). Using (3.6), we find that 

2 



A + 

3KjJ 



A + /i 



3kh 



-laij 



(3.20) 



where ij G {11,22, 12,21,23}. In particular, we recover 
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A + 2/.+ ^±^A + i? 



4/x 



1 

3k, 



A 



2A + —A - B 
2/i 



(3.21) 



for the classical acoustoelastic effect (Hughes and Kelly, 1953). Note that C12 — C21 = 1 
and that di2 — ^21 = l/(3ft), and hence in the incompressible limit di2 = c?2i- Using these 
relations, we establish that 

- ^21) = ^1 + l/(3/^)c^^ (3.22) 

Recall now that p is the mass density in the reference configuration. When expressed in 
terms of the deformed density pc = pJ~^ equation (3.22) becomes 



(3.23) 



This relation is in fact exact in accordance with (1.3) and the expressions pvf^ = JAoijij, 
i ^ j specialized accordingly with (T2 = 0. 

Finally, we take the incompressible limits of the elastic constants using the limiting 
values listed in Section [l] There, we have pvf^ — )■ 00 and pf I2 — ?• 00, unsurprisingly, be- 
cause longitudinal homogeneous plane waves may not propagate in incompressible solids. 
For the transverse principal waves travelling in the direction of tension. 



A 



pvf^ = p + I3p + - ] ei + I5p + -A + 3D ] ef 



7 



in terms of the elongation (Destrade et al. , 2010a), and 

A 



PV12 = P + 1 



A 
12p 



1 

V 



2p + 



3-^]+3D 
Ap) 



a 



l5 



(3.24) 



(3.25) 



in terms of the pre-stress. For the transverse principal waves travelling perpendicular to 
the direction of tension, 



pvli = /i + -ei + (2/i + A + 3D) el. 



A 



P^23 = /i - 3/i + - ei + 5/i + -A + 3D e^. 



7 



(3.26) 
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in terms of the elongation, see (Destrade et al. 



2010a), and, using (3.7), 



A f A\ 

8fi + 3A + — + 3D 
8/i 



l5 



cr 



1) 



(3.27) 



in terms of the pre-stress. 

For an apphcation, we use data for a sample of silicone rubber that has been subjected 
to a standard tensile test. Figure [2] displays the variation of the tensile Cauchy stress 
component di with the elongation up to a maximum stretch of about 250%, at which 
stage one end of the sample snapped out of its grip. Over that range, the TOE strain 



energy density (1.8) is not able to capture the behaviour of the sample adequately, as 



shown in Fig. |2| and is thus discarded. On the other hand, the FOE strain energy (1.9) 
gives an excellent least-squares fit, with coefficient of correlation = 0.9997. In fact, 
the fit is very good up to a much larger value of the stretch than could be expected of this 
fourth-order approximate theory. We determined the following values for the constants 
of second-, third-, and fourth-order elasticity: 



/i = 109.35 kPa, A = -454.18 kPa, D = 109.27 kPa. 



(3.28) 



Note that all three are of the same order of magnitude, as expected from the theory 
(Destrade and Ogden, 2010). We also remark that — 8/x < A < — 4/i, indicating that 



the corresponding Mooney-Rivlin solid is materially stable (Destrade and Ogden 2010 



Destrade et al. 2010b). 



Using these values, the theoretical variations of the squared wave speeds pvi2, pvii, 



pv23 are plotted versus the elongation in Fig. [3j These curves suggest that it is sufficient 
to elongate the sample by about 20% to reveal quadratic acoustoelastic effects, and thus 
to determine D experimentally. 



4 Discussion 

This study was set within the framework of the theory of small-amplitude waves prop- 
agating in a deformed elastic solid. Existing general expressions for the speeds of body 
waves in a finitely deformed isotropic elastic solid were recalled and then specialized to 
fourth-order elasticity. Specifically, the squared wave speeds were expanded in terms of 
the pre-deformation (and, equivalently, in terms of the pre-stress) in order to reveal the 
so-called acoustoelastic effect at the considered order. While the classical acoustoelastic 
effect is concerned with the linear variation of the squared wave speed with the pre-stress, 
the expansion was extended to the quadratic regime (large acoustoelastic effect). Exphcit 
expressions for the quadratic coefficients for compressible and incompressible materials 
were determined for the cases of hydrostatic and uniaxial pre-stress. These coefficients 
give access to the fourth-order elastic constants. 

In the case of an incompressible sample under uniaxial tension, for example, mea- 
surement of the variations of the speed of a single transverse wave with respect to the 
pre-stress (or pre-strain) is enough to determine the elastic constants. The wave speed 
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Figure 2: Fitting tensile data (thick curve) of the Cauchy stress ai (in Pa) versus the 
elongation ei for a sample of incompressible silicone rubber; the TOE strain-energy den- 
sity cannot capture the variations adequately; the FOE model gives an excellent fit with 
only 3 material parameters, with values fi = 109.35 kPa, A = —454.18 kPa, D = 109.27 
kPa. 



in the unstressed configuration gives the initial shear modulus n of second-order elastic- 
ity, the linear variation of the squared wave speed with the pre-stress or pre-strain gives 
the third-order Landau-coefficient A, and the quadratic variation gives the fourth-order 
constant D. 



For compressible materials, Hughes and Kelly (1953) showed that any three of the 



seven equations (2.5), (3.9), (3.14), with just the linear terms retained, or equivalently 
(2.8), (3.19), excluding either the expression for pvf2 or for pt^l^, suffice to determine the 



two Lame coefficients A and fi of second-order elasticity and the three Landau coefficients 
of third-order elasticity A, B, C. Here we have shown that any four of these seven 
equations, with the quadratic terms now included, give access also to the fourth-order 
constants E, F, G, H (again, excluding either the expression for pw^j fo^ P'^2i)- 

Although we have restricted attention largely to third- and fourth-order elasticity the 
expressions for the components of the tensor of instantaneous elastic moduli ^Aq given in 



(1.2 ) apply for an arbitrary finite deformation relative to an unstressed configuration of an 



isotropic elastic material, and the subsequent incremental response depends on the finite 
deformation and its accompanying pre-stress. If instead there is an initial stress in the 
reference configuration then the components of ,Aq are considerably more complicated. 



as detailed in Shams et al. (2011), but they do clarify how the elastic moduli depend in 



general on the initial stress. In particular, the dependence of .A^ on an initial hydrostatic 
stress T can be put in the simple form 

■Aoijki = lJ'iT'){6ik6ji + 5ii6jk) + \{T)5ij5ku (4.1) 
where ^{t) and A(r) are now stress-dependent Lame moduli, with r > (< 0) corre- 
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Figure 3: Plots of pv"^ (in MPa) versus elongation ei for a deformed incompressible block 
of silicone rubber subject to simple tension; theoretical curves for waves travelling (a) 
in the direction xi of tension {v = V12), (b) in any direction in the plane normal to the 
direction of uniaxial stress, with transverse polarization in that plane {v = V23), (c) in 
any direction in the plane normal to the direction of uniaxial stress, with polarization 
normal to that plane (f = V21); the dotted hues correspond to the linear part of the pv'^ 
versus ei curves (the classical acoustoelastic effect). 



sponding to tension (pressure), their precise functional dependence determined by the 
choice of strain-energy function. We refer to a recent paper by [Rajagopal and Sacco- 



mandi (2009) and references therein for a discussion of the stress dependence of moduli 



within the context of an implicit theory of elasticity. 
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